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ABSTRACT 



Non-Gaussian imprints on the cosmic microwave background radiation (CBR) sky are 
within the grasp of current experiments. A clear non-Gaussian signature would be point- 
like sources. We have examined the nature of possible point sources that were tentatively 
identified in a recent high frequency CBR experiment with half-degree resolution (Cheng et 
al. 1994: MSAM). The effects of local foreground sources, including cold dust clouds, ra- 
dio sources and the Sunyaev-Zeldovich effect due to foreground rich clusters are considered, 
and the effective spectral slopes of these various foreground contaminations are calculated. 
Radio source emission and the Sunyaev-Zeldovich effect are ruled out as the explanation of 
the possible MSAM sources. Modeles are examined of extremely cold dust clouds which are 
located in the solar neighborhood, the interstellar medium, the galactic halo or at cosmo- 
logical distances. We estimate the cloud mass and dust grain parameters, and in particular 
the grain size distribution, that are required in order to produce a detectable signal in an 
MSAM-type experiment. It is shown that cold dust clouds can have an important effect on 
CBR experiments only if the clouds are nearby, and located within a disc scale height of the 
solar neighborhood. Cold dust emission remains a possible source of far-infrared signal at 
the level of the detected CBR fluctuations on degree scales, but may be ruled out if the dust 
emissivity index satisfies a = 1.5 ± 0.5 . 



1 Introduction 



Valuable information about the early universe and the physical processes that generate the 
primordial fluctuations from which cosmic structure formed can be gained from measure- 
ments of cosmic microwave background radiation (CBR) temperature anisotropics (White, 
Scott & Silk 1993). An especially important issue is the Gaussian nature of the temperature 
anisotropics. Gaussianity in the linear regime is a generic conswequence of most inflationary 
theories for the origin of the fluctuations. However none of these models fare particularly 
well at accounting for the large-scale power spectrum of density fluctuations on all observed 
scales, and rather extreme solutions have been advocated (e.g. Peacock and Dodds 1994; 
Peebles 1994; Bartlett, Blanchard, Turner and Silk 1994). Hence deviations from Gaussianity 
are a possibility that can only be limited by experiment. Although significant progress has 
been made in the understanding of the CBR since the detection of fluctuations at 7° angular 
scale by the COBE satellite (Bennett et al. 1992; Smoot et al. 1992; Wright ct al. 1992), the 
Gaussianity question remain unresolved. Despite the full-sky coverage achieved by COBE, 
the combination of beam smoothing and the effects of cosmic variance (Luo 1994; Hinshaw 
ct al. 1994) preclude COBE alone from testing Gaussianity. At the same time, while the 
detection of degree-scale fluctuations intrinsic to the CBR is a remarkable achievement, the 
datasets at these intermediate angular scales are still too small to carry out Gaussianity tests, 
although tests at these angular scales would certainly be decisive once large sky coverage is 
achieved (Coulson et al. 1994). 

Despite all these difficulties, there are non-Gaussian imprints on the CBR sky that are 
within the grasp of current on-going experiments. A clear non-Gaussian signature would 
be the detection of point-like sources. In fact, two candidates for such sources may have 
been detected by the medium scale anisotropy measurement (MSAM) experiment. There 
are serious issues of data analysis that pertain to whether or not possible point sources are 
subtracted before attempting to measure temperature fluctuations: we do not address such 
issues here. Rather, we ask the question: could possible foreground sources produce point 
source-like signals in a GBR experiment at MSAM resolution and frequency? 

Various topological defects, notably soft-domain wall bubbles (Goetz & Notzold 1991; 
Turner et al. 1991), the global monopoles (Bennett & Rhie 1991) or texture (Turok & 
Spergel 1991) are capable of producing spotlike GBR anisotropics of any desired size by 
choosing appropriate model parameters. However, before relying on topological defects as 
the interpretation of candidate sources, one has to carefully filter out any foregrounds. In 
carrying out the experiments, GBR anisotropy signals have to be separated carefully from 
local foreground sub-millimeter and millimeter radiation fields. Three possible foreground 
sources are studied in this paper. These are cold dust clouds, nonthermal extragalactic radio 
sources and the Sunyaev-Zeldovich effect in foreground galaxy clusters. Multifrequency 
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measurements have previously been studied as a ttechnique for removing the foreground 
(Brandt et ah 1994), with a focus on the point-hke sources hsted above. 

The arrangement of this paper is as follows: an effective spectral index for point sources 
is introduced in section 2, and results for radio sources, dust clouds and the SZ effect are 
presented and discussed in section 3. We conclude that both radio sources and the SZ 
effect are ruled out as a possible explanation of MSAM-type sources, that is to say, several 
sigma fluctuations that are point- like at ~ 30 arc-min resolution and have a spectral energy 
distribution that, crudely at least, is indistinguishable from that of the CBR. Cold dust 
emission remains a possibility, but we find that this option also may be ruled out if a 
conventional value for the dust emissivity Q; = 1.5±0.5is adopted. 



2 Effective Spectral Slope of Point Foreground Sources 

A multi-frequency technique has been used in most CBR experiments for the purpose of 
removing foreground emission. The spectral information in each pixel of the observations 
is used as a discriminator between foreground and the true CBR signal. In the MSAM 
experiment on which we focus in this paper, four frequency channels have been used, 
u — 5.6,9.0,16.5 and 22.5cm"^. By fitting a model for cold dust emission to these four 
frequencies, we can study the spectral features of the emission, and the range of parame- 
ters which can still give rise to a spectrum that is close to that of the cosmic microwave 
background radiation. 

Let us first define the effective spectral slope of cold dust emission. Given the spectrum 
of dust emission, S^, the effective power-spectral slope a^ff is 

S. = {-r"lBM), IbM) = B^{To)^(^^) , (1) 

1^0 — 1 V i / rms 

where is a reference frequency, which we choose it to be — 5.6cm'~^, B^{To) is the 
blackbody radiation spectra, Tq — is the CBR temperature and x — hu/KTo — 

2.97^. At half-degree scales, the expected temperature anisotropy is ( ^ j ~ 1.0 x 10~^. 

V / rms 

Thus, the expected CBR flux Ib{v-,Tq) is 

lB{y,To) = 6.5 X io5(^)4e2-9W-o(g2.97./.o _ i)-2jySr-i. (2) 

By this definition, the effective spectral slope cte// = for true CBR anisotropies. 

One can determine the effective spectral slope from measuring the flux Si at four different 
frequencies i/j. Let iji — log[Si/ lB{i^i,TQ)] and Xi — log{vi/vQ), the best fit spectral slope can 
be found through minimizing x^, 

X^ = Y.{yi-A-Bxif/N, iV = 4, (3) 
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here A, B are two constants and the best fit slope B is our effective spectral slope cte//, 
which is found to be 

_ E^X^y^/N-{E^X,JN){J:,yi/N) 

We now discuss the flux S^- The flux in each pixel is sampled by a two-beam or three-beam 
chop. We will consider the situation where in one position there is a compact source of 
unknown origin. By a compact source, we mean that the angular size of the source is smaller 
than the beam width, which is 6^ = 0.425 x 28' — 12'. To simplify the problem, we consider 
a two-beam square-wave chop (in the MSAM experiment a sine-wave chop is used, and the 
position of the compact source matters if the angular size is much less than the beam width). 
In this case, the observed flux is 

S,^7rel,lB{iy,To)S + j,, (5) 

where 

\ Iq / \ Iq / rms 

is the normalized temperature anisotropy. It is a random Gaussian variable of zero mean 
and unit variance. The MSAM beam width is 9jn — 0.425 x 28' = 12 arcminutes, and thus 
the expected CBR flux at uq — 5.6cm~^ is / = 1.4Jy. 

We parametrize the observed flux from any pixel that contains a compact source by a 
parameter r, 

= 7rellB(u,To)[5 + rI(u/uo)], (7) 

where r is the ratio of the contribution from the compact sources to CBR temperature 
anisotropies and is the frequency dependence which is normalized so that I{v) = 1 

at u — vq. 

As shown in Fig.l, three possible foreground sources are studied in this paper. These are 
cold dust clouds, nonthermal radio sources and the Sunyacv-Zeldovich effect in foreground 
galaxy clusters. For each case, the flux j^, from each source is discussed in detail in the 
following. 

2.1 Radio Sources 

The flux is assumed to be 

3{v) = A{^f, (8) 

here vq — 5.6cm~^ is the reference frequency and B —1 is the spectral index of the radio 
emission from candidate sources. 
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For the radio sources as modeled above, 

/ 9d\'^ A T( I \ (\(\r.^(^ \B-i -2.97 -i^/ 2.97-i^-lx2 

^ = i:^' ^ ' -(e ^0 ). (9) 

Here 9d is the angular size of the radio source. 



2.2 The Sunyaev-Zeldovich Effect 

The typical angular scale of a foreground rich cluster is of order several arcminutes, which 
lies in the range of the MSAM sources. The scattering of microwave photons by hot electrons 
in the intracluster gas will make a rich cluster a powerful source of submillimeter radiation. 
The flux density is given by (Sunyaev & Zcldovicli 1980) 



■ e^ + l ; 
X - 4 



MTo), (10) 



1)' 



where x — hu/kTo and y — J{kTe/'meC^)aTnedl. 

For the SZ effect, the ratio r and frequency dependence /(— ) is 



2.97v 



.Ujn/ [OJ-/J-)rms T^Q Vl^Q/e —1 

2.3 Dust Emission 

Let us consider a single dust cloud which is located a distance D away from our location. 
The flux from the dust cloud is (Greenberg 1968; Kwan & Xie 1992) 

3. = 1 1 nd{rd,T)B,{T)e,7iA7irle-^'^drddT, (12) 

where By{T) = 4.63 x 10^i/^(exp ^z/ — l)~^Jy5'r~^ is the brightness function of a black-body 
with temperature T and v is in units of 100 GHz. Here, naird, T) is the distribution of dust 
grains as a function of grain size ra and grain temperature T, — ^{^2^)°' emissivity 
of the dust grains and t^, is the optical depth of the dust emission. Since reabsorption of 
the dust emission is expected to be negligible, the emission should be optically thin so that 
= 0. When the cloud consists of dust paricles of a single radius, and therefore of a single 
temperature 7^, the flux from the dust cloud is: 

^ -^'^^l^^B^i^d), (13) 
here, Nd is the total number of dust particles in the cloud. 



5 



For dust emission, in addition to the parameter r which denotes the ratio of dust emission 
flux to the expected rms CBR flux, we have two additional parameters: the dust temperature 
Trf and the emissivity index a. r and /(^) are given by 



3 Results and Discussions 

At the four frequencies Ui = 5.6, 9.0, 16.5, 22.5cm~^, -i = 1, 4 , which are used in the MSAM 
experiment, the observed flux is parametrized as 



The physical interpretation of r is straightforward: r/(l + r) is the percentage of the total 



take into account the fact that 5 is a random variable. In our calculations, we generate 100 
Monto-Carlo realizations of 6 each time when computing aeff, and the mean and variance 
of for the different foregrounds are shown in Fig. 2 (a-f). 

3.1 The Sunyaev-Zeldovich Effect 

In Fig. 2a, the effective spectral slope for the SZ effect is shown. Again, we find that at 
the la level, r greater than 0.30 is ruled out based on spectral analysis. We conclude that 
no more than 23% of the signal comes from the SZ effect of a foreground rich cluster if the 
spectral slope lies within ±0.3 of the CBR spectral index. 

3.2 Radio Sources 

In Fig. 2b, the effective spectral slope for B = —1.0 radio sources are shown for different 
ratios r. Again at the la level, r greater than 0.40 is ruled out based on spectral analysis. 
We conclude that no more than 28% of the signal comes from radio emission if the spectral 
slope lies within ±0.3 of the CBR spectral index. 

3.3 Cold Dust Emission 

Estimation of the dust emission is more complicated than the calculation of the radio source 
contribution and the SZ effect. We need to specify the dust temperature T^, the dust 
emissivity a and the geometry of the dust grains. The dust geometry may be dominated by 
whiskers or fractals: we will address this point later. In Figs. 2c, 2d, 2e, 2f, the effective 




(14) 



S,^ = 7c91Ib{i^„ To)[6 + rI{u,/uo)],z = 1, 4. 



(15) 



flux contributed by foreground sources at 5.6cm ^ To determine cte// from Eq.(^), one has to 
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spectral slopes for Td = AK cold dust with dust emissivity a = 1.5, 1.0, 0.5, 0.0 are shown. 
The spectral slope of dust emission can be close to that of the CBR if the dust emissivity 
is small [a = 0.5,0.0), as shown in Figs. 2e, 2f. Dust emission can be responsible for the 
putative MSAM point sources if the dust emissivity is small. However, if we adopt for the 
physical value of the dust emissivity a = 1.5 ± 0.5, then no more than 50% of the total flux 
can be due to cold dust. 

How much dust can give rise to a flux that could cause problems for medium-scale 
CBR experiments? For the time being, let us assume that the MSAM 'sources' are due 
to foreground dust emission, and we will estimate the required dust mass. To determine the 
cloud mass, we choose the mass density of dust grains to be pd = 3g/cm^ (Hildebrand 1983). 
From Eqs. (13), the mass of the dust cloud is 



M = %rlNd = 8.58 x io28(^)(^£1^)2(^)(^)1-«(^)2,-3-.[^m.. _ ^jg^^^^ 
3 IJy 6 e 1/im Ipc 

(16) 

Two 'sources' are reported at z/ = 5.6cm~^ by the MSAM experiment (Cheng et al. 1994). 
The flux from one source is 3.7±0.9 Jy, the other is 2.9±0.7 Jy. Both sources are unresolved 
at 6py/HM = 28'. Let us assume that the angular size of the dust cloud is Od = 14' in order 
to minimize the beam-smoothing effect and dust emissivity a ~ 1.5. The estimated mass is 



Ml = (1.2 ± 0.3) X 102^(^^)(— ) Vams (17) 

1/im Ipc 



for the 3.7 ± 0.9 Jy source and 

Ms = (0.9 ± 0.3) X 102^(^^)(— )Vams (18) 

1/im Ipc 

for the 2.9 ± 0.7 Jy source. 

The dust clouds may be located in the solar neighborhood {D < 100 pc), the interstel- 
lar medium ( 100 pc < D < Ikpc), the dark halo (5kpc < D < 50kpc) or even be at 
cosmological distance {D ~ 3000 Mpc). Let us look into all these possibility in detail. 

3.3.1 Dust in the Solar Neighborhood 

Here, the solar neighborhood is loosely deflned as the region centered around the sun and 
within the disc scale height {D < 100 pc). In the inner solar neighbourhood, D < Ipc, the 
dominant heating source is the sun. The dust mass budget is fairly small (M < lO^^g). 
However, the typical dust temperature is too hot to account for the MSAM sources. The 
dust temperature is 
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where e ~ 0.1 for dielectric dust grains, a ^ 1.5 is the gain emissivity index and a is the grain 
radius. For typical dust grains of radius = 0.5/im, the grain temperature is = 50K for 
D = 0.1 pc. 

In the far solar neighborhood, D ~ 100 pc. One expects the heat source to be diffuse star 
hght. Dust can be very cold if the radius of the particles is large (Rowan- Robinson 1991) or 
if they are fractals or needles (Wright 1993). Let us first consider fractals and needles. The 
physics of far infrared emission from fractals or needles differs from those of the spherical 
dust we considered in the previous section. Since fractals and needles have cooling times 
which are shorter than the average arrival times between incoming photons, the emission 
from fractals and needles are not determined by their equilibrium temperature but rather 
by their enthalpy U{T), mass spectrum and absorption efficiency (Wright 1993). 

The problem with explaining point sources with emission from fractal dust or needles 
is that the spectra of the IR emission differs from the CBR spectrum considerably due to 
the transient heating by single photon absorption events. The heat capacity of the dust is 
very small at low temperatures (T < lOK), and absorption of a single optical photon will 
increase the temperature to a very high degree unless the fractal or needle building block is 
large enough. The average energy of a photon in the radiation field given by Eq. (2) is 

E = SkTs^ = 2.3ev (20) 
C(3) ^ ' 

for Tg = lO^i^. The enthalpy is 

rr ^rr.^ 4.15 X lO-^^T^'^ 

^^^^^ = (1 + 6.51 X 10-3r+ 1.5 X 10-^r^ + 8.3 X lO-^T-^f ^'/'''"^ ^^'^ 

for graphite (Guhathakutra & Draine 1989). The temperature increase of a building block 
consisting of N atoms when absorbing one photon of energy E is 

Thus, > 10^*^ to avoid excessive spectral deviation from CBR, or, the radius of building 
block must be greater than ~ O.lfim, which is about the size of typical dust grains that 
absorb in the optical band. For fractals of dimension D, the emissivity is enhanced by 
{L/ri,Y~^ , where L is the size of the fractal grain. To reach an equilibrium temperature of 
several Kelvin, L must be much greater than ri,. Thus, even for fractal grains the size is 
large. 

For clouds of large radius (r^ ~ SO/^m) grains to produce MSAM-like sources, the dust 
mass is 

Ma - 1Mq(^—Y. (23) 
°400pc^ ^ ' 
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Taking a typical gas-to-dust ratio rj = 160, the mass of the gas cloud is M = AOMq for 
D fti 50pc. The angular size of the clouds located at D = 50 pc is 10' for / = 0.15 pc. 
Although it is very unlikely that dust can be cooled down to a few Kelvin in the solar 
neighbourd, this window remains open for explaining the MSAM sources. 

3.3.2 Clouds in the ISM 

In this case, the clouds are concentrated in the galactic plane. There will be no high galactic 
latitude dust emission from these dust clouds and hence there will be no effects on the MSAM 
experiment, which samples high galactic latitudes, even if they exist in the ISM. 

3.3.3 Dust in Cold Halo Gas Clouds 

It has been suggested that cold molecular gas clouds might be the dominant form of galactic 
dark matter, distributed either in a large disk (Pfennigcr, Combes and Martinet 1994) or 
flattened halo (Gerhard & Silk 1994). The latter model leads to the natural possibility of cold 
dust clouds far from the galactic plane. The dust temperature at the periphery of the galaxy 
will be as cool as a few Kelvin even if the conventional, physically acceptable dust emissivity 
is adopted. One initial worry is that dust emission from such clouds at high galactic lattitude 
would have dramatic effects on the experiments that are searching for microwave background 
fluctuations at millimeter and submillimetcr wavelength. We will show here that for the cloud 
models given by Gerhard & Silk (1994), there is virtually no effect on the GBR from dust 
emission. In the model where the clouds are in near-hydrostatic equilibrium, for a cloud of 
density lO'^cm^^ and temperature 100 K to satisfy the Jeans criterion requires the size and 
mass of the cloud to be L = Aj/2 = 0.64pc and M = 85A^23^^0- Here A^23 is the hydrogen 
column density in units of 10^^ cm~^. Demanding the cloud-cloud collision time scale to be 
at least a Hubble time constrains the column density to be Nh > 8 x 10^^ cm^. Thus, the 
mass of the clouds is hmited to be M < IOMq. Even if dust dominates in these clouds, 
the flux at 5.6cm~^ will be Si, ~ ImJy. The angular size of the dust region is 9 a ~ 10". 
After smoothing over the beam, the observed flux will be around 1//Jy, several orders of 
magnitude smaller than the flux from CBR temperature fluctuations. 

3.3.4 Primeval Dust 

The last possibility of high latitude dust emission that we consider arises from primeval 
galaxies at high redshift (Bond et al 1986). The angular size of the dust envelope around a 
primeval galaxy at redshift z is given by: 



l{l + z) 
Dh 



1 



)■ 



(24) 
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Here, / is the linear size of the primeval galaxy, which is around 100 kpc. The typical angular 
size of the possible point sources produced by a primeval galaxy is: 

/ z 

9 = 1.14/i(— — )(-)arcminute (25) 

The temperature of the CBR scales linearly with redshift z, Tc — 2.73(1 + z)K. At the 
periphery of the primeval galaxy with D=50 kpc, the dust temperature is — 2A{^)K for 
grain size = 0.5//m, which is hot in the usual sense. The observed flux density will be: 

S,^s = 3.45Jy(M/7.85 x lQ^MQ){—^f{hQkpc/DHf{l + zf+^ (26) 

"fwhm 

To explain the MS AM experiment in the framework of primeval dust, the dust mass should 
be: 

Md ~ lO^^M©. (27) 

This is about three orders of magnitude larger than the typical dust mass in observed IR 
galaxies. The bolometric luminosity of the MSAM sources if they are located at cosmological 
distances is L ~ lO^^L©, much brighter than the brightest-known IRAS galaxy F10214+4724. 
It is improbable that the MSAM sources are primeval galaxies. 

In conclusion, we have examined the nature of the putative MSAM sources and shown 
that foreground radio sources or and the SZ effect in foreground rich clusters are ruled out as 
dominant contributions. Unless the dust emissivity is unphysically small, cold dust emission 
can also be ruled out. 

This work is supported in part by the NSF. 
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FIGURE CAPTION 

Fig. 1: The spectral density of various sources at submillimeter range. The solid line is the 
CBR anisotropics; The dotted line is SZ flux from a rich cluster; The short dashed line is 
for radio sources ( spectra index -1); The long dashed lines are for 4K cold dust emissions 
with different dust emissivity a = 1.5, 1.0, 0.5, 0.0, from the top to the bottom. 
Fig. 2a: Best fit spectral slope for SZ effect. 
Fig. 2b: Best fit spectral slope for radio sources. 

Fig. 2c: Best fit spectral slope for 4K dust emission with dust emissivity a — 1.5. 
Fig.2d: a = 1.0. 
Fig.2e: a = 0.5. 
Fig.2f: a = 0.0. 
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